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ABSTRACT 

We present follow-up observations of two galaxy clusters detected blindly via the Sunyaev- 
Zel'dovich (SZ) effect and released in the Planck Early Release Compact Source Cata- 
logue. We use the Arcminute Microkelvin Imager, a dual-array 14-18 GHz radio interfer- 
ometer. After radio source subtraction, we find a SZ decrement of integrated flux density 
— 1.08 ±0.10 mJy toward PLCKESZG12 1.1 1+57.01, and improve the position measurement 
of the cluster, finding the centre to be RA 12 h 59 m 36.4 s , Dec +60°04' 46"8, to an accuracy of 
20". The region of PLCKESZG1 15.71+17.52 contains strong extended emission, so we are 
unable to confirm the presence of this cluster via the SZ effect. 

Key words: cosmology: observations - radiation mechanisms: non-thermal - cosmic mi- 
crowave background - galaxies: clusters - Sunyaev-Zel'dovich - galaxies: clusters: individual 
(PLCKESZ G 1 1 5 .7 1 + 1 7.52, PLCKESZ G 1 2 1 . 1 1 +57 .0 1 ) 



1 INTRODUCTION 

The Planck satellite is an instrument designed to measure the full 
Cosmic Microwave Background (CMB) sky to the highest sensi- 
tivity yet achieved. Part of its remit is a blind survey of galaxy 
clusters via their Sunyaev-Zel'dovich effects (SZ; Sunyaev & Zel- 
dovich 1970); this occurs when CMB photons are inverse-Compton 
scattered by hot gas in galaxy clusters. The resulting shift in the 
blackbody spectrum of the CMB produces a decrement at radio fre- 
quencies below the null of 217 GHz and a positive signal above this 
frequency (see e.g. Birkinshaw (1999) and Carlstrom et al. (2002) 
for reviews). 

The Arcminute Microkelvin Imager (AMI) is a dual-array ra- 
dio interferometer operating at 14-18 GHz; part of its scientific 
programme is making observations of hundreds of known galaxy 
clusters to measure their gas masses and structures. Along with its 
high sensitivity, its resolution of ~ 2'.5 allows the measurement of 
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the extent and shape of galaxy clusters above a redshift of ~ 0. 1 
(below which they become resolved out). 

In 2011, the Planck Early Release Compact Source Cat- 
alogue (ERCSC; Planck Collaboration 2011b) was made pub- 
licly available; it includes a sample of 189 galaxy clusters de- 
tected via their SZ effects at signal-to-noise (S/N) > 6 (Planck 
Collaboration 2011a). AMI was used to follow up the two pre- 
viously unconfirmed clusters lying in the northern sky: PLCK- 
ESZG115.71+17.52 (S/N= 6.78) and PLCKESZG121. 11+57.01 
(S/N= 6.66), AMI observing details of which can be found in 
Tab. 1. 



2 AMI SZ OBSERVATIONS 

AMI comprises two aperture synthesis arrays located at the Mullard 
Radio Astronomy Observatory, Lord's Bridge, Cambridge, UK. 
The AMI Small Array (SA) consists of ten 3.7-m diameter 
equatorially-mounted dishes with a baseline range of ~ 5-20 m, 
while the AMI Large Array (LA) has eight 12.8-m-diameter dishes 
with a baseline range of ~ 20-100 m. Both arrays observe I+Q in 
the band 14-18 GHz, each with system temperatures of about 25 K. 
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Name 


RA 


Dec 


SA observing 


SA map noise / 


LA observing 


LA 19-pt map 


LA 61-pt map 




(J2000) 


(J2000) 


time / hrs 


,LtTy beam - 1 


time / hrs 


noise/ /iJy beam" 1 


noise/ [iJy beam - 1 


PLCKESZG115.71+17.52 


22 26 24.9 


+78 18 15.8 


50 


78 


77 


55 


200 


PLCKESZG121. 11+57.01 


12 59 23.8 


+60 05 24.8 


40 


81 


67 


65 


250 



Table 1. Statistics of the AMI observations of the two candidates. 



Channel 


v/GHz 


5 ^2«6 /Jy 




FWHM LA /' 


3 


13.9 


3.74 


1.89 


6.08 


4 


14.6 


3.60 


1.78 


5.89 


5 


15.3 


3.47 


1.68 


5.70 


6 


16.1 


3.35 


1.60 


5.53 


7 


16.9 


3.24 


1.52 


5.39 


8 


17.6 


3.14 


1.45 


5.25 



Table 2. Assumed I+Q flux densities of 3C286 and 3C48 over the 
commonly-used AMI band, and the LA primary beam full-width-half- 
maxima (FWHM) for each channel. 



The backends are analogue Fourier transform spectrometers, 
from which the complex signals in each of eight channels of 750- 
MHz bandwidth are synthesised, and the signals in the synthesised 
channels are correlated at the ~ 10 per cent level. In practice, the 
two lowest-frequency channels are generally not used due to a poor 
correlator response in this frequency range and interference from 
geostationary satellites. Further telescope details are given in AMI 
Consortium: Zwart et al. (2008). 

The AMI data reduction was performed using our in-house re- 
duction software REDUCE. This is used to apply path-compensator 
and path-delay corrections, to flag interference, shadowing and 
hardware errors, to apply phase and amplitude calibrations, make 
a correction for the changing intervening air mass and to Fourier 
transform the correlator data to synthesise the frequency channels, 
before outputting to disk in wv-FITS format suitable for imaging in 
the Astronomical Image Processing System (AIPS; Astronomical 
Image Processing Software 2007). Flux calibration was performed 
using short observations of 3C48 and 3C286 near the beginning and 
end of each run, with assumed I+Q flux densities for these sources 
in the AMI channels consistent with Baars et al. (1977) (see Tab. 2). 
As Baars et al. measure I and AMI measures I+Q, these flux den- 
sities include corrections for the polarization of the sources. The 
phase calibrators, J0017+8135, J1302+5748 and J2133+8239, were 
chosen from the Jodrell Bank VLA Survey (JVAS; Patnaik et al. 
1992) on the basis of proximity « 5°) and flux density (^ 0.7 Jy at 
15 GHz). After phase calibration, the phase of AMI over one hour 
is generally stable to 5° for channels 4-7, and to 10° for channels 3 
and 8. The overall consistency of the flux-density scale is estimated 
to be better than five per cent. 

To analyse the SA cluster observations we use a Bayesian 
analysis package, McADAM (Marshall et al. 2003 and Feroz et al. 
2009b). This uses MULTlNEST (Feroz & Hobson 2008 and Feroz 
et al. 2009a), an implementation of nested sampling (Skilling 
2004), to efficiently explore the multidimensional parameter space 
and to calculate the Bayesian evidence. This analysis has been ap- 
plied to pointed observations of known clusters (e.g. AMI Consor- 
tium: Rodriguez-Gonzalvez et al. 2010; AMI Consortium: Zwart 
et al. 2010), and also to detect previously unknown clusters (AMI 
Consortium: Shimwell et al. 2010). We use it here to model the ef- 
fect of radio point sources on the S A data, using the LA data as pri- 
ors. To transform the SA data into a format suitable for McADAM, 
the Mv-data for all good observations were concatenated together 
to make a wv-FITS file for each channel. These were then trans- 
formed into lists of visibilities for the purpose of generating a co- 
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Figure 1. LA raster signal-to-noise maps of the two candidates, with 
sources detected at 4o~ and above marked; their properties can be found 
in Tables 3 and 4. Contours are S/N= 1, and follow levels +3, +4, +5, +6, 
+8, +10 and +20. The small ellipse in the box in the lower-left shows the 
average synthesised beam. 



variance matrix for the data. The covariance matrix describes the 
terms which contribute to the data but are not part of the model, i.e. 
the CMB, thermal noise from the telescope and confusion noise 
from unresolved point sources. The latter were modelled using the 
Tenth Cambridge Radio Survey (10C: AMI Consortium: Franzen 
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et al. 2010; AMI Consortium: Davies et al. 2010), integrating the 
confusion power from zero up to the LA source detection limit. 

The data were also binned into cells of width 40 A, determined 
in Marshall (2003) as a suitable binning scale to reduce the size 
of the data to a manageable level without adversely affecting the 
resulting inference of model parameters. The binned data files were 
then analysed by McADAM using the method described in AMI 
Consortium: Zwart et al. (2010). Maps were made using IMAGR in 
MPS from each channel of the SA and LA observations; however 
we here present only the combined-channel maps of the SA and LA 
observations. The task imean was used on the LA individual maps 
to attach the map noise to the map header, and f latn was then used 
to stitch the maps together, with a primary beam correction applied 
using parameters shown in Tab. 2. 

Source-finding was carried out on the LA data using software 
modified from that used in the 9C survey (Waldram et al. 2003). 
Spectral indices were fitted using LA maps of all six channels, as- 
suming source fluxes follow a power-law relation of S °= v~ a for 
the AMI frequencies. The properties of point sources detected at 
4 o~ and above by the LA are used as priors when modelling the 
SA data in McADAM. The position prior used in McADAM is a 
delta function since the resolution of the LA is around three times 
that of the SA; two exceptions are discussed in Section 3.1. LA flux 
densities generated the peak of Gaussian priors on the SA source 
flux densities, with as set to the calibration error of five per cent 
of the source flux density added in quadrature with the local map 
noise. Spectral index (a) priors were also set as Gaussians with a 
equal to the error on the spectral index fit. These errors tend to be 
small (~ 0.4) for bright sources and large (~ 2) for faint sources. 
Fig. 1 shows the LA maps for the candidates with 4o"la ra dio point 
sources marked; the positions and properties of the sources are 
given in Tables 3 and 4. 



3 RESULTS AND DISCUSSION 

After the source models were fitted to the SA data in McADAM, the 
sources were subtracted from the data using the modelled source 
fluxes and spectral indices. The data were then mapped again in 
AIPS to generate source-subtracted maps. The parameters used for 
source-subtraction are the last two columns of Tables 3 and 4. 

3.1 PLCKESZG1 15.71+17.52 

The raw map for this candidate, Fig. 2a, shows a very complicated 
region with both positive and negative features. The sources ap- 
pear to consist of two main complexes: one to the south-east of the 
pointing centre, consisting of sources L, I, J, N, K, M and O. To the 
west and north is a separate complex, consisting of sources D, P, 
C, B, G and the extended source A. Given the extent of this source, 
and the uncertainty on whether source B was also extended, these 
were modelled with Gaussian position priors with a = 5". This is 
useful when dealing with extended emission as the location and 
magnitude of the peak of such emission may differ significantly 
when measured on different scales; it has previously been used 
to account for small amounts of extended emission in AMI Con- 
sortium: Shimwell et al. (2011) and AMI Consortium: Rodnguez- 
Gonzalvez et al. (2011). 

After source-subtraction (see Fig. 2b), extended positive resid- 
uals are visible directly to the north of the pointing centre, and to the 
south of source A. These areas of emission have no obvious coun- 
terparts on the LA map (Fig. la) and have SA peak flux densities 



~ 0.4 mjy. However the complexes have integrated flux densities 
of ~ 0.7 mjy each. This seems to indicate that the emission is real, 
extended, and positive. 

A difficulty arises when using an interferometer to measure 
structure of unknown extent, as it is unclear how much will be 
resolved out due to the lack of wv-coverage on large scales. This 
partially negates the natural filtering advantage of using an interfer- 
ometer for SZ detection, discussed in Birkinshaw (1999). When de- 
convolving the dirty map with the synthesised beam in the CLEAN 
process (Clark 1980), incomplete interferometric measurement of 
extended structure can lead to artifacts, including negative signal. 
It is unclear whether the negative signal 7' south-west of the point- 
ing centre, in the source-subtracted map, is purely artifactual from 
poor CLEANing of the extended positive emission, or also includes 
a component from SZ emission. Publicly available ROSAT data 
(Voges et al. 1999) show elongated X-ray emission of ~ 5'extent 
at the Planck ESZ position. If this object is a galaxy cluster then 
the SZ and X-ray emission are likely not to be strongly offset, so 
the south-west negative feature (and other negative features on the 
map) are probably artifactual. 

As to the source of the positive extended emission seen in this 
region, we note that Galactic observations with the SA (see e.g. 
AMI Consortium: Hurley-Walker et al. 2009) tend to pick up more 
extended emission closer to the Galactic plane. This is likely to be 
synchrotron associated with regions which also have bright thermal 
dust emission. The latter was also noted by Planck Collaboration 
(201 la) as a problem when measuring the SZ effect toward this ob- 
ject. There is also the possibility that the two source complexes are 
extended radio doubles from jet-producing galaxies, which would 
be more likely to appear in a galaxy cluster than in the field. 

3.2 PLCKESZG121.11+57.01 

It is clear even before source subtraction that there is isolated excess 
negative signal in the centre of the SA map, Fig. 3a. After source- 
subtraction the initially elliptical appearance is revealed as a close 
to circularly-symmetric SZ decrement, Fig. 3b. The 6o" positive 
residual between sources P and N corresponds well to an overden- 
sity of 3o" features nearby in the LA map, and is likely to be ex- 
tended emission resolved-out by the LA. The feature is more than 
one synthesised beam away from the cluster and appears alone, so 
should not affect measurement of the cluster decrement. 

From the source-subtracted map, we measure the integrated 
flux density of the central negative signal as — 1.08 ± 0.10 mjy. 
We estimate the S/N of our measurement as 13, by dividing the 
absolute value of the integrated flux density by the map noise 
of 0.08 lmjy. Thus we confirm the Planck blind SZ detection of 
PLCKESZG121. 11+57.01 as a galaxy cluster. We hope to inves- 
tigate the full comparison of AMI and Planck data via common 
cluster models in future publications. 

We also measure a position offset from the ERCSC po- 
sition, finding A RA = +94 ± 13" and A Dec = -38 ± 19", re- 
sulting in an improved cluster position of RA 12 h 59 m 36.4 s , 
Dec +60°04' 46f'8. Wen et al. (2009) identify a cluster of redshift 
0.33 at RA 12 h 59 m 33.4 s , Dec +60° 04' 09" in the Sloan Digi- 
tal Sky Survey DR6 Galaxy Clusters Catalog. Given that the sep- 
aration between these two measurements is merely 44", and the 
redshift is typical of clusters detected by Planck (Planck Collabo- 
ration 201 la), it is likely that the two are the same object. Publicly- 
available ROSAT data (Voges et al. 1999) show only a few X-ray 
photons, consistent with noise, so the cluster appears to have a low 
X-ray luminosity. 
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ID 


RA 


Dec 


S LA /mJy 


a LA 


SMcADAM/mJy 


•^McAdam 


A 


22 23 41.0 


+78 18 35.1 


6.73±0.83 


0.55 ± 0.21 


8.71 ±0.09 


0.79 ±0.11 


B 




+78 19 49.4 


0.50 ±0.23 


0.48 ± 1.55 


0.50 ±0.09 


0.41 ± 1.27 


C 


22 23 53.6 


+78 21 29.1 


3.64 ±0.61 


0.43 ± 0.45 


4.32 ±0.08 


0.68 ±0.19 


D 


22 25 29.4 


+78 25 25.3 


2.53 ±0.5 1 


0.89 ± 0.53 


2.47 ±0.08 


1. 14 ±0.28 


t 


22 23 41.6 


+78 10 51. 1 


2.02 ±0.40 


1.14 ± 1.58 


2. 18 ±0. 12 


2. 13 ±0.60 


F 


22 21 12.6 


+78 13 08.6 


1 .75 ±0.61 


— 0.99 ± 1.77 


1.03±0.35 


— 1.16± 1.68 


Q 


22 24 35 3 


+78 1 8 05 5 


1.67 ±0.41 


2.13 ± 0.61 


1 80 + 07 


1.25 ±0.33 


H 


22 21 09.2 


+78 21 40.8 


1.37±0.47 


0.37 ±1.73 


1.63 ±0.23 


1.56±1.13 


I 


22 26 41.7 


+78 17 39.7 


1.14 ±0.34 


0.74 ±0.70 


1.16±0.09 


1.67±0.39 


J 


22 26 36.9 


+78 16 53.3 


1.12±0.34 


0.28 ±0.72 


1.35 ±0.09 


0.97±0.39 


K 


22 26 43.9 


+78 14 08.0 


0.82±0.29 


0.63 ±1.01 


0.74 ±0.07 


0.04 ±0.74 


L 


22 26 28.4 


+78 19 55.8 


0.77 ±0.28 


-0.23 ± 1.01 


0.09 ±0.05 


-0.01 ±0.98 


M 


22 27 34.7 


+78 15 00.6 


0.72 ±0.28 


0.25 ± 1.04 


0.91 ±0.07 


0.13±0.81 


N 


22 26 34.7 


+78 15 24.9 


0.31 ±0.18 


1.15 ± 1.56 


0.41 ±0.06 


0.90±1.38 





22 27 59.3 


+78 16 31.5 


0.28 ±0.1 8 


0.90 ± 1.70 


0.16 ±0.06 


0.73 ±1.62 


P 


22 24 50.3 


+78 22 57.7 


2.06 ±0.42 


0.30±0.73 






Q 


22 26 49.9 


+78 24 28.7 


0.53 ±0.27 


0.70±1.37 






R 


22 23 56.9 


+78 15 12.8 


0.47 ±0.23 


0.23 ± 1.66 






S 


22 30 06.4 


+78 14 34.4 


0.40 ±0.40 


0.40± 1.63 







Table 3. Sources found in the LA map of PLCKESZ Gl 15.71+17.52, in order of modelling prior and LA flux density. Sources A and B were modelled with a 
Gaussian position prior of a = 5" to make allowance for the extended emission detected in this area. Sources P, Q, R, and S were directly subtracted from the 
SA data to reduce the computational resources needed to model the sources in the centre. 



ID 


RA 


Dec 


S LA /mJy 


«LA 


SMcADAjn'mly 


"mcAdam 


A 


12 57 58.3 


+59 53 32.9 


3.97±0.79 


-1.08 ± 1.35 


2.91 ±0.37 


-0.49±1.15 


B 


13 00 15.4 


+60 12 01.3 


2.66 ±0.55 


0.20 ±0.95 


2.55±0.13 


-0.12±0.55 


C 


13 00 40.7 


+59 55 50.6 


2.47 ±0.61 


0.22 ± 1.46 


3.26 ±0.24 


0.52±0.84 


D 


13 01 05.9 


+60 10 49.4 


1.76±0.58 


-1.04 ± 1.63 


1.13±0.23 


-0.02 ±1.47 


E 


12 59 03.0 


+59 50 38.7 


1.71 ±0.57 


— 1.48± 1.69 


1.26 ±0.28 


-1.68±1.55 


F 


12 58 36.2 


+60 00 04.8 


0.97±0.34 


0.64±1.39 


0.69 ±0.10 


0.65 ±1.26 


G 


13 00 36.8 


+60 07 58.1 


0.85 ±0.35 


-0.99 ±1.53 


0.51 ±0.12 


-0.92±1.39 


H 


12 5941.6 


+60 03 25.2 


0.70 ±0.27 


— 0. 17 ± 1.13 


0.74±0.12 


1.27 ±0.75 


I 


12 58 45.3 


+60 01 49.0 


0.68 ±0.27 


-1.54 ± 1.37 


0.58 ±0.09 


-1.92± 1.22 


J 


12 59 29.6 


+60 02 31.8 


0.64 ±0.26 


-0.77 ± 1.12 


0.63 ±0.09 


-0.25 ±0.95 


K 


12 58 42.0 


+6007 39.5 


0.47 ±0.23 


1.34 ± 1.56 


0.57 ±0.08 


1.67 ± 1.27 


L 


12 58 39.6 


+60 05 45.0 


0.46 ±0.23 


0.46 ± 1.66 


0.32±0.08 


0.66 ±1.48 


M 


12 59 52.9 


+60 01 33.2 


0.41 ±0.21 


1.19± 1.53 


0.52 ±0.08 


0.89± 1.38 


N 


12 59 42.9 


+60 11 29.1 


0.38±0.21 


0.68 ±1.75 


0.48 ±0.09 


0.09 ± 1.47 


O 


12 59 39.6 


+60 04 08.9 


0.34 ±0.19 


— 1.22± 1.64 


0.34±0.09 


-0.76 ±1.47 


P 


12 59 03.1 


+60 08 43.8 


0.29±0.18 


-0.79±1.71 


0.19±0.07 


-0.51 ±1.64 


Q 


12 59 51.2 


+60 08 30.3 


0.28±0.18 


-0.77 ±1.57 


0.44 ±0.07 


0.41 ±1.33 


R 


13 00 03.2 


+59 48 03.9 


4.68±1.35 


-0.47 ± 1.77 






S 


13 01 46.8 


+60 07 49.5 


3.66±1.04 


-0.74 ± 1.72 







Table 4. Sources found in the LA map of PLCKESZ G121. 11+57.01, in order of modelling prior and LA flux density. Sources R and S were directly subtracted 
from the SA data to reduce the computational resources needed to model the sources in the centre. 



4 CONCLUSIONS 

We confirm the Planck blind SZ detection of PLCK- 
ESZG121. 11+57.01 as a galaxy cluster at 13rj, with a 15GHz 
integrated flux density of —1.08 ± O.lOmJy. Our observa- 
tions improve the measurement of the cluster position to 
RA 12 h 59 m 36.4 s , Dec +60° 04' 46"8. We cannot positively 
identify PLCKESZG115.71+17.52 as a galaxy cluster via SZ, as 
the data are consistent with both the presence of an SZ decrement 
and extended positive emission, or simply the latter. However we 
note that radio source environment may indicate the presence of 
jets from active galaxies, often found in galaxy clusters. 
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Figure 2. PLCKESZG1 15.71+17.52 SA map before and after source subtraction; contours are every O" = 78/iJy, and follow levels —20, —10, —8, —6, —5, 
—4, —3, —2, +2, +3, +4, +5, +6, +8, +10 and +20. Sources marked with a '+' were modelled in McADAM, while those marked with a ' x ' were directly 
subtracted. The sources marked with a A were given a Gaussian position priors of width 5" due to their slightly extended structure. The parameters of the 
labelled sources can be found in Tab. 3. The small ellipse in the box in the lower-left shows the synthesised beam. 




RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 

(a) Before source subtraction. (b) After source subtraction. 

Figure 3. PLCKESZG121. 11+57.01 SA map before and after source subtraction; contours are every O" = 8 1 /X Jy, and follow the same levels as in Fig.2; 
annotations and markings are also the same. Sources E, R, and S are outside the field-of-view shown here, and the latter two were directly subtracted from the 
SA data. The parameters of the subtracted sources can be found in Tab. 4. 
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